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Abstract. We investigate two models for the following setup: We consider 
a stochastic process X S C[0, 1] whose distribution belongs to a parametric 
family indexed by i? 6 8 C 1. In case 1? = 0, X is a generalized Pareto 
process. Based on n independent copies X^ , . . . , X (™) of X, we establish 
local asymptotic normality (LAN) of the point process of exceedances among 
XM, . . . , X ( n ) above an increasing threshold line in each model. 

The corresponding central sequences provide asymptotically optimal se- 
quences of tests for testing Ho : i3 = against a sequence of alternatives 
H n : -& = # n converging to zero as n increases. In one model, with an un- 
derlying exponential family, the central sequence is provided by the number 
of exceedances only, whereas in the other one the exceedances themselves con- 
tribute, too. However it turns out that, in both cases, the test statistics also 
depend on some additional and usually unknown model parameters. 

We, therefore, consider an omnibus test statistic sequence as well and com- 
pute its asymptotic relative efficiency with respect to the optimal test sequence. 



1. Introduction 

In the recent three decades, the focus of univariate extreme value theory shifted 
from the investigation of maxima (minima) in a sample to the investigation of 
exceedances above a high threshold. This approach towards large observations 
eased accessing the field of extreme value theory and became a crucial tool for 
various applied disciplines, such as building dykes. 



Since the publications of the articles by Balkema and de Haan (1974) and 



Pickands ( 1975 ) it is known that exceedances above a high threshold can reasonably 
be modeled only by (univariate) generalized Pareto distributions (GPD), resulting 
in the peaks-over-threshold approach (POT). 

Due to practical necessity, the focus of extreme value theory moved in recent 
years to multivariate observations as well. Accordingly, the investigation of multi- 
variate exceedances enforced the definition and investigation of multivariate GPD. 
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This investigation is still lively continuing as even the definition of multivariate 



GPD is under debate; see, for instance, Falk et al. (2010 Chapter 5) 



As already mentioned by de Haan and Ferreira (2006 p. 293): Infinite- dimen 



sional extreme value theory is not just a theoretical extension of multivariate ex- 
treme value theory to a more abstract context. It serves to solve concrete problems as 
well. Such concrete problems are, e.g., observing dykes and tides along their whole 
width and not only at a finite set of observation points. There is, consequently, the 
need for a POT approach for functional data and for generalized Pareto processes 
as well. Again, the data exceeding some kind of a high threshold are modeled by 



a functional counterpart of a GPD; see Aulbach et al. (201261. The current paper 



deals with optimal tests that check for particular models whether those exceedances 
do, in fact, arise from such a corresponding process. 



1.1. Previous work. Following Buishand et al. (2008), a standard generalized 



Pareto process, i.e., a generalized Pareto process with ultimately uniform tails in 
the margins, is defined as follows. For convenience, we use bold font such as V for 
stochastic processes and default font such as / for non stochastic functions. All 
operations on functions such as / < are meant pointwise. 

Definition 1.1. Let U be an on (0,1) uniformly distributed random variable (rv) 
and let Z — {Zt)te[o.i} G C[0> 1] be a stochastic process on the interval [0, 1] having 
continuous sample paths. We require that U and Z are independent and choose an 
arbitrary constant M < 0. Then 



(1) 



V := (V t ) te[0A] 



max I , M 

Zt 



te[o,i] 



defines a standard generalized Pareto process (GPP) if < Z t < m, E(Z t ) = 1, 
t £ [0,1], hold for some constant m > 1. A stochastic process Z € C[0, 1] with 
these two properties will be called a generator. 

The constant M is incorporated in the above definition to ensure that Vt > — oo 
for each t € [0,1]. Note that the finite-dimensional marginal distributions of V 
provide multivariate GPD with ultimately uniform tails; see, e.g., Aulbach et al.| 



(2012a 



The process V is characterized by the fact that its functional distribution func- 
tion (df ) is given by 



P{V <f) = l-E\ sup (|/(t)|Z t ) 
\te[o,i] 



/G£-[0,1], H/IU < aro, 



for some x > 0. We set E~[0, 1] := {/ € E[0, 1] : / < 0} where E[0, 1] denotes 
the set of all bounded functions / : [0, 1] — > K that have only a finite number of 
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discontinuities. Then 



l/ll 



E 



sup (\f(t)\Z t ) 

,i£[Q,l] 



/e£[o,i], 



defines a D- 



on -E[0, 1] with generator Z; see Aulbach et al. (2012c). This 



representation of the df of V in terms of a D-norm is in complete analogy with the 



multivariate case of a GPD. We refer to Falk et al. (2010 Section 5.1 



It was established by Aulbach et al. (2012c I that a stochastic process X £ 



C[0, 1] is in the functional domain of attraction of a max-stable process (MSP) 
£, denoted by X € if and only if this is true for the univariate margins 

together with convergence of the corresponding copula processes. A stochastic 
process U = (U t )te[n.i] G C[0, 1] is called copula process if U t is for each t G [0, 1] 
uniformly distributed on (0, 1). 

For each standard GPP there is a corresponding standard MSP, i.e., a stochastic 
process t] = {vt)te[o.i] € C[0, 1] such that 

(2) P(77</)=exp(-||/|y, feE-[0,l]. 

Note that this implies P(r?t < x) = exp(x), x < 0, t G [0,1]. On the other hand, the 
df of each max-stable process r\ having standard negative exponential margins has 
a representation as in Equation ([2]); we refer to Aulbach et al. (2012c) for details. 

1.2. Overview of the current paper. We replace the rv U in Equation ([!]) by 
a rv W > which is independent of Z, too. However, the distribution of W is 
different from the uniform one and, thus, the process 

W 



(3) 



X 



( x t)te[o, 



-,M 
Zt //te[o,i] 



is no longer a standard GPP. 

This gives rise to the following problem: Based on the exceedances in a sample 
of n independent copies X^\ . . . , X^ 71 ' of X above a high threshold line, how close 
can the df of W get to that of U with the difference still being detected? As we 
consider exceedances above a high threshold, only the lower end of the df of W 
matters. In other words, the problem suggests itself to define parametric models 
{H$ : $ £ 0} for the df H$ of W, such that we can derive optimal tests detecting 
the deviation of the distribution of the upper tail of X from that of V, i.e., the 
deviation of t? from zero. This is the content of the present paper, which is organized 
as follows. 

In Section [2] we require that the df H$ of W has a density h$ near zero, which 
satisfies for some 5 > the expansion 

(4) h § (u) = l + $u s + o(u s ) as u I 
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with some parameter i)e 6, where zero is an inner point of 9 C 1. The standard 
exponential df, for instance, satisfies this condition with 3 = 1 and $ = — 1. The 
null-hypothesis i? = is meant to be the uniform distribution on (0, 1). 

In Section [3] we assume that the distribution of W belongs to an exponential 
family given by the probability densities on the interval 

(5) h#(u) = C{d) exp(i?T(u)), < u < 1, & € R. 

In both models wc establish local asymptotic normality (LAN) of the point 
process of exceedances among . . . , X^- n ' above an increasing threshold line. 

The results, which are stated in Theorem |2.2| and Theorem |3.2| provide in each 
model the corresponding central sequence and, thus, optimal tests for testing d = 
against a sequence of alternatives i? n converging to zero as the sample size increases. 
It turns out that the particular values of the exceedances contribute to the central 
sequence only in model Q, whereas in the exponential model ([5| the number of 
exceedances alone yields the central sequence. The fact that just the number of 
realizations in shrinking sets provides the central sequence was characterized for 



truncated processes in quite a general framework in Falk ( 1998) and Falk and Liese 



(1998) 



As the central sequences and, thus, the asymptotically optimal tests also depend 
on further parameters of the generator process Z , which might be unknown in 
practice, we consider an omnibus test for testing -d = as well. We compute its 
asymptotic relative efficiency (ARE) with respect to the optimal test in each model. 
While ARE is positive in model Q, it is zero in model 

2. Testing in ^-neighborhoods of a standard GPP 

This section deals with optimal tests in the model introduced in Q . We assume 
that the df of the rv W > in ([3| belongs to a parametric family {H^ : § £ 0} of 
distributions, where 9 is an open subset of K containing 0. By H we denote the 
uniform distribution on the interval (0,1). In addition to Q, it is required that 
there is some uq € (0, 1) such that the density h#(u) of H${u) exists for u e [0, uq], 
i)g9, and satisfies for some 5 £ (0, 1] the expansion 

(6) h § (u) = 1 + -du s +r#(u), ue[0,u o ], 

where r#(0) = 0, i? € G, and 



(7) sup 

0<|1?|<£0 



du 5 



0(1) 



as u l for some eo > 0. Obviously, Q is equivalent with r#(u) = o(i9u s ) as u \. 0, 
uniformly for |$| < Eq. Since we have h Q = 1 and r = 0, Q and ^ imply in 
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particular the representation 



M") = ho(u) ll + O[(H (u)) 5 



i.e., the lower tail of H§ is in a 5-neighborhood of H$; see Falk et al. (2010 Sec- 
tion 2.2). 

Moreover, we assume 



(8) A := El inf Z t > 0. 

V ' \te[o.i] 7 

As inf t6 [ 01 ] Z t > 0, this condition is equivalent with the assumption that inf tg r j] Z t 
is not the constant function zero. ^ is, for instance, satisfied if Z = 2U, where 
U = (Ui) te r X ] € C[0, 1] is a copula process such that — U e 2? (77), 77 e C[0, 1] being 
a standard MSP. This is implied by the fact that in this case P (inf te [o,i] Ut > 0J = 
1. 

Note that ^ and Holder's inequality also give 

B:=E[ inf Z}+ 5 ) > 0. 
V*e[o,i] y 

2.1. Local asymptotic normality. In order to derive asymptotically optimal 
tests in this model, we first establish local asymptotic normality (LAN) of the 
point process of exceedances 

N n , c (B) := £ e ( (i)/c) (fl H [0, 1]), B£ B, 

where X^ l \ i < n, are independent copies of X in ([3| and c < 0. B denotes the 
(T-field of Borel sets of R and £ K is the point measure with mass one at x. Note that 

(9) sup — < u <^=> X >cu, ue [0, 1], 
te[o,i] c 

i.e., the random point measure N niC actually represents those observations among 
X^\ . . . , which exceed the constant threshold function c. 



Denote those observations among sup te [ ,i] (x^/cj with sup tG [ 0il ] (x^/cj 
1, i < n, by Yi, . . . , in the order of their outcome. Then we have 

N n , c (B)= J2 £ Y k ( B )> Bel. 

k<r(n) 



< 



By Theorem 1.4.1 in Reiss (1993) we may assume without loss of generality that 
Yi, Y2, . . . are independent copies of a rv Y with df 

my<u) = p ;%I c ^ o<«<i, 

P&(X > c) 

under parameter and that they are independent of the total number r(n), which 
is binomial B (n, P$ (X > c))-distributed. 
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In the next lemma we provide the density /^ jC of sup 4g [ -n (Xt/c) and, thus, the 
density of Y under i9, which is f$ :C /P$(X > c). By P*Z we denote the distribution 
of a rv Z. 

Lemma 2.1. Suppose that the distribution of the rv W in ([3| belongs to the family 
{Htf : & £ 0}. Then there is some cq < such that the density fy tC ( u )> w*^ 1 respect 
to the Lebesgue measure, of the rv sup 0<t<1 (X t /c) exists for § € © } c € [co,0), 
u G [0, 1], and it is given by 



UA U ) = \ C \ / zhs(\c\zu) P* inf Z t ) (dz). 
Jo V *G[Q.i] / 

Furthermore there exists eo > smc/i i/ia£ 

/#, c (u) = |c| A + |c| 1+5 B« J + o (tf |c| 1+a ) 

uniformly for < eo awe? it e [0, 1] as c f 0; note i/iaf foA u ) = l c l 

Proof. Let m be given as in Definition 1 1 . 1 1 and uq, £q, S be given as in Equation 
Then we obtain for Co > max {M, — u /m}, $ £ <d, c £ [co,0) and u € [0,1] by 
conditioning on inf te r 0;1 i Z t = z and Fubini's theorem 

P#(X > cu) = P$ (max [ --^r,M ] > cu, i € [0, 1] 



Z t 

Po(W< \c\uZ u t£ [0,1]) 



Ptf [W < \c\u inf Z* 
V te[o,i] 

^ Pd(W < |c| zw) ^P * ^ inf^ Z t J (dz) 



fl^flclzu) P* inf ZA (dz) 
te[o,i] 



|c| / z ft,,j(|c| zx) ( P* inf Z t ] (dz)dx. 



This representation implies that sup t6 r 0il i (X t /c) has the density 
f§A u ) = \ c \ I zhi)(\c\zu) ( P* inf Z t I (dz) 

Jo V ^[O. 1 ] / 

= \c\ J z(l + tf(|c| zit) 5 + r^(|c| zu)) ^P * ^ inf^ Z^ (dz) 

- |c|£ U f .i z ') + *" c|1+< " l£ ((4i f ., z ') 1+i ) 

= |c| A + \c\ 1+S u s B + o(ti |c| 1+5 ) 
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as c | 0, uniformly for < e an d u € [0, 1]. (Note that r#(\c\ zu) = o yd |c| j 
uniformly for \d\ < £q, u £ [0, 1] and z £ [0, m] as c t 0.) As Ho(u) — u, we have 
h (u) =1, a£ [0,1], which completes the proof. □ 

If Co < c < 0, we have Pq(Y < u) = u, i.e., the distribution C$ tC {Y) of Y under 
-d is then dominated by Cq iC {Y). Thus the distribution £$(N n>c ) of N n ^ c under $ 
is, in this case, dominated by £o(N n c ) - see, e.g., Theorem 3.1.2 in Reiss (1993) - 
and we obtain from Reiss (1993 Example 3.1.2) the £o(-^n,c)-density of C^(N n>c ) 

dU(N n , c ). . _ U, c ( yi ) P (X > c) 



d£o{N n 



>)= n 



t\ foMP*(X>c) 



P*{X>c) 
Po(X > c) 



M([0,1]) 



I -PAX >c) 
1 - P (X > c) 



n- M ([0,l]) 



where fi = Y^i3i s Vi , < y±, . . . , y M r[o,i]) < 1 an d A*([0, 1]) < n. The loglikelihood 
ratio is, consequently, 

d&>(N n 



(10) 



= E ** 

i<r(n) 



dCo(JV„,c) (JVn '' 

/^(Fi) P (X > c) 



+ (n - r(n)) lo. 



/ ,c(F) Ptf(X > c 

1 - P4X > c 



-(n) log 



^(X >c) 
P (X > c) 



l-P (X>c) 

We let in the sequel c = c n depend on the sample size n with c n j" and, equally, 
$ = with i?„ — > as n — > oo. Precisely, we put with arbitrary £ € R 



(ii) 



:=0«(£) 



25 



1/2- 



The following result finally proves the desired LAN property of A„, c ; it is a crucial 
tool for deriving asymptotically optimal tests in the subsequent subsection. 



1+25 



Theorem 2.2. Suppose that c n f 0, n\c„ 
for i? n as in (11) f/ie expansion 

Ln^iK I 0) = {1 + s)A1/2 (Znl + Z n2 ) - 2A{2S + 1) 



The 



OFo(l) 



obtain 



£ 2 P 2 



with 



2A(2S + 1) ' A(2S + 1) 
r(n) - n |c„| A 



(n|c n |A)V2 



^Z5„ #(0,1), 
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and 



Z n 2 



1 + 5 

r(n)V2 



fc<r(n) 



>£>„ iV 0, 



J 2 



2(5 + 1 



6emg independent. 



Proof. First we compile several facts that will be used in the proof. From Lemma 
12.11 we obtain 

(Fact 1) P Q (X > c n ) = P(V > c n ) = \c n \A 

and, thus, a suitable version of the central limit theorem implies 

r(n) - n |c„| A 



(Fact 2) 



(n\c n \Ay/2 



■ N(0,1). 



Moreover, we conclude from Lemma 2.1 for |i9| < Eq 

(Fact 3) P4X > c n ) - P (X > = \c n \ 1+S Y~§^ + ( tf l c «| 1+5 ) 

and, thus, 



(Fact 4) 



P?„(X>c„)-P (X>c n ) 



1 



1 



Po(X>c„) (n|c n |)V2(i + ^ VMcnl) 1 / 2 

Hence, Taylor expansion log(l + e) = e — e 2 /2 + o(e 2 ) as e — > implies 



r(n) log 



^„(-y > C„) 

P (* > C„) 



+ (n — r(n)) log 



l-P? n (X>c n ) 
1 - P (X > Cn) 



PyJX > C n ) - Pp(X > C n ) _ 1 fP,„(X • r„) - ,n,(_Y • r„ 



Po(X > c„ 



P^(X>c n )-P (X>c„) 



P (X > Cn) 



Po{X>c n ) 

Pp{X > C n ) - Py n (X > C n ) 

1 - P (X > c n ) 

+ O (\P# n (X > Cn) - P (X > C TI 



+ (n - r(n)) 



, sPtiJX > Cn) - Po{X > Cn) . ,,P (X > C n ) - P § JX > C„) 

T ( n ) 5 , v ^ I \ + (n - r(n))- 



P (* > c„) 
r(n) /P?„(X > c„) - P (X > c„) 



1 - Po(X > c n ) 



2 V 



Po{X > c 



OP (1) 



r(n) 



P^„(X>c„)-P (X>c„) 



Po(^ > c„) 



P (r(n))O 



(n|c n |)3/ 2 



nP (X > c n )0 



(n|c„|)3/ 2 



0, 
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where ~ denotes asymptotic equivalence. The preceding convergence to zero follows 
from the condition n |c„| 1+2<5 — > n _>oo 00 an d the equivalence 

(n - r{n)) \P* n {X > c n ) - P (X > c n )\ 2 ~ O (|c n |) 

From the law of large numbers we obtain 

( P# n (X > Cn) - P (X > c 

r(n) 



£(r(n)) 
71 IcJ A 



P (X > On) 

1 £P ( 1 



(71^1)1/2(1 + 5)^ V(^KI) 1 / 2 

1 £P / 1 



(n|c n |)V2(l + ( 5)^ \(n\c n \y/i 
t 2 B 2 



™ (l + tJ) 2 ^ 

Moreover, 

P# n (X > c„) - P (X > c„) > c„) - P# n (X > c n ) 

r v n ) p c v ^ + (n - r(n > 

P (A > c„) l-Po(X>c„) 

P 1?n (X>c n )-P (^>c„) 



P (X>c n )(l-P (X>c n )) iT{n) - nPoiX>C ^ 
{n\c n \A)^( 1 iL> 1 



(n\c n \)^(l + S)A(l + o(l)) V>KI) 1/2 
r(n) - nP (X > c ra ) 
(nK|A)i/2 

^nl+Op (l). 



(1 + <5)AV2' 
Altogether we have shown so far that 

f P* n (X>Cn) \ , , , ^1-P^(X>C„) 



p (^>c„)y v K " °U-p«(i>c„) 
cp _ t 2 B 2 

Z n \ - Wf. ; .s 2 . + Op (l). 



(1 + 5)AV2 ni 2(1 + 5) 2 A 
Next we show that 
v-* , ( h n .cAYk) P (X>c n ) \ _ jB j 2 B 2 5 2 

t fr/° S V/o,c n (^) P# n {X>c n )) ~ Ai/ 2 (l + S) Zn2 2A(2S + l)(l + 5) 2+0p ^ ) - 

We have by Lemma pO] 

/0,c„(l1fe) {n\c n \) L / 2 A 

where ro(Yfc, c„) = o ^(n |c„|) 1 ^ 2 ^j uniformly for fe and n with 
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f 1 ( U n ,Cn(t) _ . _ g B g\ h, Cn {t) 

o V /o,c n (i) (n|c„|)V2^ ; P Q (X>c n ) 

P#JX>c n )-P (X>c n ) £ B [ x tSdt 



P (X>c n ) {n\c n \y/ 2 A J 

P# n (X > c n ) - P (X > c n ) £ B 



and 



Po(X>c n ) (n\c n \y/ 2 (l + 6)A 

Var (r (Fi, n , c„)) < P (rg(H, c„)) = o (l/(n |c„|)) . 



Using again the Taylor expansion log(l+e) = e — e 2 /2 + 0(e 3 ) as e — > 0, we deduce 



fe<r(n) 

— r(n) lo; 



P^(X>c»)-Po(X> CTt ) 



_ t \( PvAX > c n ) - P (X > c n ) _ 1 / P^(X> CTI )-P (X> Cn ) ^ J " 
TW 1 V Po(X>c„) "" 2 V P (X>c n ) 

+ op (l) 

^ 2 v y" i T(n) ^ i of m 

2n|c n |A 2 ^ fe n|c„|2(l + <5) 2 A 2 PoW 



fc<r(n) 

r(n) 1 / 2 1 

(n| C „i)V2 T ( n) i/2^: fe ^ v ifc 1 + * 

r(n) £ 2 P 2 1 v2 * , r{n) £ 2 P 2 

k<r(n) 



2n|c„| A 2 T („) 2-, r * + „ | c „| 2(1 + <5) 2 ^ 2 + P ° l J 



& i ^ / 4 i \ es 2 s 2 

2^ \ Y k TT7 9479ATTT7TXAW + ° PoUj 



AV2 T ( n )V2 ^f,V 1 + <V 2^(25 + 1)(1 + 5) 2 
-+d N[-- 



2A(2S + 1)(1 + S) 2 ' A(25 + 1)(1 + 
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by the law ol large numbers and the central limit theorem. This completes the 
proof. □ 

2.2. Testing z9 = against <d = #„. Denote by u a = $ _1 (1 - a) the (1 - a)- 
quantilc of the standard normal df. By the Neyman-Pearson lemma and Theo- 
rem |2.2| the test statistic 

/ (25 + l) 1 / 2 \ 

fl ( N n,c n ) = l(«„,oo) [ j— j (Znl + Z n2 )\ 

defines an asymptotically optimal level-a test, based on N niCn , for Ho : $ = 

/ 1+2(5 \ 

against d n — $ n (Q — £./ [ n \ c n\ ) with £ > 0. As ip 1 (N ntCn ) does not 
depend on £, this test is asymptotically optimal, uniformly in £ > 0. 

The corresponding uniformly asymptotically optimal test for Hq against 
with £ < is 

(N \ i ((26 + 1)^ , 

V?2 W V ri,e,J — L(-oo,-u a ) I \J^S 11 ' 

The asymptotic power functions of these tests are provided by Theorem |2.2| as 
well. By LeCam's third lemma we obtain that under "d n = 

Ln,c n (#n I 0) = {1 + s)A1/2 (Znl + Z n2 ) 2A{2S + 1) + Op« n (1) 



V \2A{25 + 1)' A{25 + l) 



Znl + -^ri2 —>Da„ N 



with 

gg(l+j) (1 + ^) 2 ' 
^ 1 /2( 2( 5 + 1) ' 25 + 1 
The asymptotic power functions of ifi are, consequently, given by 

(12) Hm E Km ( W ( AT„, C J) =!-*(«*- ^ 1/2( ^ 1)1/2 ) • 

A disadvantage of the optimal test statistics ipi(N n ,c n ) is the fact that they 
require explicit knowledge of the constants A and 8. To overcome this disadvantage, 
we consider in the following an alternative test. 

Recall that the observations Y%, Y 2 , . . . are independent and, under ■& — 0, uni- 
formly on (0, 1) distributed rv if the threshold c is close to zero. Conditional on the 
assumption that there is at least one exceedance, i.e., conditionally on r(n) > 0, 
the test statistic 

r(n) 



v ; fe=i 



is under Ho exactly N(0, l)-distributed. By $ we denote the standard normal df. 
This test statistic is analogous to that in Falk and Michel ( 2009 ) for testing for a 
multivariate generalized Pareto distribution. 
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The next result provides the asymptotic distribution of T nCn under the alterna- 
tive d n — asn->oo. 

Proposition 2.3. Under the assumptions of Theorem \2.S\ we have 

T n ,c n ->D.„ N U-^p J I($W)V(X) dx, lj . 

Proof. First we compute the asymptotic mean and variance of <I> -1 (Y) under i? n 
and c„ for n — > oo. From Lemma |2.1| we obtain that the density of Y under d n is 
for < u < 1 and c n > Cq given by 

U n ,c n («) 



P# n (X > c n ) 

K\ 



zh$ rl (\c n \uz) I P* inf Z t I (dz) 



P#„{X >c n ) J " V 
From Fubini's theorem and the substitution u t— > we, therefore, obtain 



1 

-1/ 



= / $ (u)po niCn (u) du 

— - / z ^~ 1 (u)h i) (\c n \uz)du (P* inf Z t ) (dz) 
Cn) Jo Jo V °^ t ^ 1 / 



P?„ (X > c r , 

\Cr, 



tn />oo 



z \ xh$ n (\c n \ <fr(x)z)ip(x) dx P* inf Z t I (dz) 



P# n (X >Cn) J J_ 00 V 0<t<l 

where (^(x) = $'(a;) = (27r)~ 1 / 2 exp(— x 2 /2), x € R, is the density of the standard 
normal df 

From condition ^ we obtain the expansion 



m />oo 



z / x/i^^dcnj ^(x)z)ip(x) dx [P* inf Z t )(dz) 



f zj x (l + n (|c| $(x)z) 5 + r*„ (|c„| p(a:) (p * Q haf i Z t j (dz) 

= $ n \c n \B [ x(^(x)) S cp(x)dx + oU n \c n \ S ). 
From [Fact ll and lFact 31 we obtain 

P$ n (X > Cn) = Po(X > Cn) + \c n \ 1+S + O U n | C „| 1 

l + o V 
= \c n \ A + |c„| 1+5 + o (d n \c n \ 1+t ' 

and, thus, 



1 1+5 



P„„, c „ {$-\Y)) 



TESTING FOR A GENERALIZED PARETO PROCESS 13 



0„ \c n f B x (<Z>(x)) S <p(x) dx + o (d n \c n \ 5 ) 
A+\c n \ 5 {f- s B + o(i) n \c n \ 5 ) 

j^tt.B x (<P(x)) 5 V {x) dx + o \c n f) 
A + \c n \ S ^- s B + o(# n \c n \ S ) 

1 x ($(x)) da; + o 



Equally, we obtain 

= pjx> Cn) r z sz x2kM Hx)zMx) dx ( p * & Zt ) {dz) 

I I p rn /-ex) 

— / z/ x 2 (l + & n (\c n \*(x)z) s + rt n (\c n \*(x)z))x 

n) Jo J-oo V ' 



1 



x 95(2;) dx ( P * ^ ) 



and, thus, the asymptotic variance of is under i9„ and c„ for n — > oo 

equivalent to 1. 

Finally we have the expansion 

E# n , Cn (T(n)) = nP# n (X > Cn) 

= n \c n \ A + (n \ Cn \fl 2 ^- 5 + o (n \c n \ 1+S ) 

= n\c n \A(l + o(l)). 

Now we can compute the asymptotic distribution of T n>Cn under i9 n . We have 

r(n) 
r(n) 



where the first term is by a suitable version of the central limit theorem asymptot- 
ically standard normal distributed, and 

rW^X.c, (S^QO) ~ ^„,c n (r(n) 1 / 2 ) E^ n>Cn [*~\Y)) 

~ (n M A)V" ^ ^ |)i/2 - J ^ x {${x)) s <p{x) dx 
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R f°° 

"^SiTa/ x{${x)) s <p(x)dx, 

J — oo 

which completes the proof. □ 
From Proposition |2.3| we obtain that 

( N n,c n ) ■= l(u<»,oo) (Tn,c n ) , V?2 ( N n,c n ) ■= l(-oo,-u a ) (+n,c„) 

are one-sided tests for testing $ > and i? < 0, respectively, against 0. Their 
asymptotic power functions are given by 
0(0 := lim (p* (TV )) 

(!3) / B f°° x , \ 

= i - * («a - iei ^iTa y x($(x))vw <k) , e e e. 

The asymptotic relative efficiency of if* (iV nCn ) with respect to ^ (JV„ Cn ) is, by 



2 

5, 



= (2(5 + 1) / x(<f>(x)y<p(x)dx 



(12 1 and ( |13[ ), given by the ratio 

|£| Bf^xWxWvWdx/AVf 2 
(|£|5/ (AV2 (2( 5 + i)i/2)) 2 
which is independent of £. 

Denote by fc„ := minjfc e N : ^„(f) (-^fe.cj) > ^„(?) (v< (^VcJ)} the least 
sample size, for which (p* (-/Vfc„, C(c „) is, at # n (£); a t least as good as ipi (N n , Cn ). The 
relative efficiency of <^>* (Nk„, Ch ) with respect to y>j (N n . Cn ) is then defined as n/k n . 
From ( 12 ) and ( 13 ) we obtain that 

(14) lim "I^Cl = ( 2S + V*) d * 

n ~^°° k n \c kn \ \J-oo 



see Section 10.2 in Pfanzagl (1994) for the underlying reasoning. This explains the 
significance of the asymptotic relative efficiency defined above. 

3. Testing in an exponential family model 

In this section we assume that the distribution of W belongs to an exponential 
family given by the probability densities on the interval [0, 1] 

h#(u) = C{d) exp(i?T(u)), < u < 1, i? € E, 

where T : [0, 1] — > E is a bounded Borel-measurable function satisfying 



limT(u) =: C G E, 



and C($) is defined by 



:= — , i9 g 

/ exp(i9T(u)) du 
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Figure 1. Asymptotic relative efficiency as in (14 1 



Remark 3.1. From the arguments in the proof of Lemma \2.1\ we obtain that the rv 
sup 0<t<1 (X t /c) has for c < close to zero and each on [0, 1] the Lebesgue- 

density 



U,c( u ) = l c l / ^MM zu) [P* inf Z t {dz) 



0<t<l 



< u < 1. 



In what follows we put with arbitrary £ € 
K :=1?n(f) 



where we require that C ^ /„ T(u) dtt. 

Theorem 3.2. Suppose that \c n \ — > 0, n \c n \ — > 00 as n — » 00. XTien we obtain for 



the loglikelihood ratio in ( 10 ) the expansion 



L n ,c n (# n I 0) =£Z nl - y +0 Fo (l). 



Proof. Again we compile several facts first. 

(Fact 5) C(0„) = l-0 n / T(u)du + o(tf„), n€ 

Jo 
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This follows from the expansion exp(x) = 1 + x + o(x) as x — > 0: 

1 



J exp($T(u)) dtt 
1 

~ /„ 1 + t>»T(u) + o(i? n ) du 
1 

1 + 0„ Jq 1 T(u) du + o(0 n ) 

= l-0„ / T(ti)du+o(tf„). 



(Fact 6) P*„ (X > cO - P pT > c„) 



^ n |c„|A C- / T(u)d U + o(tf n |c n |) 





n / \ \ n 

This can be seen as follows. From Remark |3. II and [ Fact 51 we obtain 

P# n {X > c n )-P {X > Cn) 

h n ,c n (u)du- I f . Cn (u)du 

m c 1 



JO 



z I hfl n {\c n \ zu) — ldu IP * inf Z t I (dz) 



0<t<l 



z[ C(d n )exp('& n T(\c n \zu))-ldu (p* inf Z t ) (dz) 

i; JO V 0<f<l / 



P * inf Z t (dz) 

0<t<l 



»i 

|c»K n (C- / T(x)dx)A + o(\c n \0 n ), 







which is lFact 6l 

|Fact 6| together with |Fact 1| yields 

fF , 7] P^(X>c n )-P (X>c n ) 1 g 

1 J Po(X>c„) (n|c„|)i/2^i/2 1 ^( n | Cn |)i/2 

Repeating the arguments in the proof of Theorem |2.2| one shows that 
T(n) log + (n - r(nj) log ' 



p (x>c„); v v 77 °u-^><=») 



^ v Jo.c„(Yk) -Ptf„(A>c„w 
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,, T(n)-nP (X> C „) £ 2 

_e (n|c„|A)v 2 y + ^o(i) 

= £,z nl - — + o Po (i). 

It remains to show that 

'/*„,*, (n) ^o(^>c„ 

fc<-r(n) 

Repeating the arguments in the proof of |Fact 6] we obtain 

fn n ,c n { u ) ~ fo,c n (u) 
fo,c n (u) 

= jf Q z (MM zu) -1)(p* Q in£ i Z t ) (dz) 

= l/ z(C(^exp(0 n T(|c„|zu))-l) fp^inf^ (dz) 

= jJ Q zUnC-$nJ^ T(x)dx + o(tf n )^J (p * (jinf^t) (ds) 
= O(0 n ) 

uniformly for w € [0, 1] and n G N. The expansion log(l + e )=e- e 2 /2 + O (e 2 ) 
for e — > together with |Fact 7 thus, yields, 

V . f fs n .cjY k ) P (X>c n ) 
Ji, ° g \focAY k ) P^(X>c n ) 

k<r(n) 

U n , Cn (Y k ) - f 0iCn (Y k ) 
fo,cjY k ) 



EM 1 " hcn{Yk ) 



fe<r(ri) 

-r(n) log ( 1 + 

V Po(^ > c„) 

v / /tf„, e „(>fe) - k.cSYk) i / / tf „,c„(n) - /o. c „(y fc ) N 2 ' 
k t {n) \ w*w _2 v " /o,c„(^) 

, P^„(X > cn) - Pp(X > Cn) r(n) ( P# n (X > c n ) - P (X > c n ) 
T[1l) Po(X > c„) ~ 2 \ P (X > c n ) 

+ o Po 1 ] 



\n |c„| 1/2 / 
Note that 



E P<) {h n , Cn (Y))= f U niCn (u) = Po n (X>c n ) 
Jo 



and 

f , Cn (u) = \c n \A = P (X>c n ). 
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We, thus, obtain 

'h n ,c n ( Y k) Po(X>c 7 



H log 

fc<r(n) 



fo,o n (Y k ) P# n (X > 



Cr 



^ \Cn\A 2^1 ~ \Cn\A 

k<r(n) 1 k<r(n) V 1 

r(n) f P„ n (X > c n ) - P (X > c n )\ 2 { Q f 1 ^ 



2 

I„ - II„ + III„ + P , 



p (x>c„) ; °V»| c „i 1/2 , 
i \ 



' \n\c n \ 1/2 J 
From [Fact 71 we obtain 

n^i / 1 g / 1 \\ 2 £ 

t j 2 1^1)1/2 ^1/2 + °| v(ri | Cn | ) l/2jJ ^ 2 - 

Next we show that I n = op D (l). This assertion follows, if we show that 



I I ! . IV, I I ' . f V > / ■ I ^ 

(17) E B 



h n , Cn {Y k ) - p^{x > c n ) y \ = r i 

IcJA / / \n\cr. 



By elementary arguments we obtain 

„ ( { U n ,c n (Y k ) - P# n (X > Cn^ 2% 



C„| A 



A 2 



I z J h# n (\c n \zY)-h#J\c n \zu) (p* o ^Z t }{dz)\ 
E Po ((J z J exp(d n T(\c n \ zY)) - exp(tf„T(|c n | zu)) du 



\2 / / rm i>1 



P* ) (dz 



which is (17). 
Finally we have 



E 



Pa 



U n ,cSY)-\Cn\A 



2 s - 



c n \ A 

^ z(h^(\c n \zY)-l) (p*^Ztj (dz?) j 

I z(C(d n ) exp(tf„T(|c„| zY)) - 1) (p * q M ^ Z^j (dz)) 
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^E Po ^l ^|l-tf B | Tui'Mlu + <>(,)„) ) (\.-^I„T[\c„\ ;Y) +„(,)„))- I 



P * inf Z t (dz) 



■■$1 (c- jf'rwduj +o(| c „| 2 ^). 

The law of large numbers implies 

e 

TT — ^ , — — 

s - s -n 'n^oo 2 



in probability, and, hence, ( 16 1 yields 

III n -II„ = o Po (l). 



We, thus, have established (15 1, which completes the proof of Theorem 3.2 □ 
The test statistic 

<t>i(N niCn ) := l( tlQj00 )(Z ri i) 

defines by the Neyman-Pearson lemma an asymptotically optimal level-a test, based 
on N ntCn , for the null-hypothesis ■& — against the sequence of alternatives d n — 
MO = £/ ((n^DV^Va (C- /*T(u)du)) with £ >0. As <MA„, C J does not 
depend on £, this test is asymptotically optimal uniformly in £ > 0. 

The corresponding uniformly optimal test for $ = against with $ < is 

4>2(N n , Cn ) :— l(_oo,-« Q )(-^nl)- 

From LeCam's third lemma we obtain that under i?„ = i? n (£) 

in,c„(^n|0)=^„l-! r + Op»(l)^„ ^fCe 2 



2 

with 

Znl ^D* n N(£,l). 

The asymptotic power functions of <pi, i — 1,2, are, consequently, given by 
lim E P {4>i(N nt0n )) = 1 - $K - i = 1,2. 

Next we compute the performance of the statistic T„ jC = r(n) -1 / 2 J2k=i 
for the testing problem # = against $ n (£). 

Lemma 3.3. We /laue 



As T„ )Cn —>£>,, c A(0, 1), Lemma 3.3 implies that the test statistic T n , Cn is not 
capable to detect the alternative i? n = #„(£)• 
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Proof. We have 

= / $~ 1 (M)p 1?n:Cn (w) du 



z I ha (\c n \uz)du P* inf Z t (dz) 
Jo V °<t< 1 / 

z / xh$ n (|c„| <&(x)z)ip(x) dx IP* inf Z t )(dz). 



P# n (X > c n ) J J V o^ 4 ^ 1 

= \Cn\ 

Pa n (X>c n )J ~ J_ ao ~"" VnKI ~ ni -v~>~>r\~'— \^ ' <t<i 
|Fact 7| implies 



m />oo 

2 



sft-^^dcnl ^(x)z)ip(x) dx [P* inf Z t j(dz) 



m />oo 

z y a;C(i?„)exp(^„r(|c„| §(x)z)) ip(x) dx ^P^inf^ j u/:) 

m roo / rl 

z xll-tin T(u) du + o(d n ) ) (1 + d n C + o{d n ))<p(x) dx 



P * inf Z t (dz) 
o<*<i M y 



JTi /"CO 



z y x I 1 + I C - J T(u) duj + o(0„) 1 p(a:) ^ ^P * Q inf i Z t 

= o(0 B ). 
We have, moreover, 

= P /yL ^ z r x 2 h# n (K\Hx)zMx)dx IP* inf Z t j (,/:). 

P?„(A>C„)J J.^ V 0<i<l 

where 

m poo / \ 

z J x 2 htf n (\c n \ 3>(x)z)cp(x) dx fp* inf Z t \ (dz) 

= A + i?„ ( c ~ J T ( u ) du ) A + °(^«)- 
From lFact ll and lFact 61 we obtain 

/I I \ 1/2 / / I | \ 1/2N 

and, thus, 
and 

Var^^ ($- 1 (Y))=E K , Cn ((^(Y)) 2 ) - ^„ >Cn (^(Y))* 
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= l + o0?„). □ 
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